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ABSTRACT 

Accretion from a disk onto a collapsed, relativistic star - a neutron star or 
black hole - is the mechanism widely believed to be responsible for the emission 
from compact X-ray binaries. Because of the extreme spatial resolution required, 
it is not yet possible to directly observe the evolution or dynamics of the inner 
parts of the accretion disk where general relativistic effects are dominant. Here, 
we use the bright X-ray emission from a superburst on the surface of the neutron 
star 4U 1820-30 as a spotlight to illuminate the disk surface. The X-rays cause 
iron atoms in the disk to fluoresce, allowing a determination of the ionization 
state, covering factor and inner radius of the disk over the course of the burst. The 
time-resolved spectral fitting shows that the inner region of the disk is disrupted 
by the burst, possibly being heated into a thicker, more tenuous flow, before 
recovering its previous form in ~ 1000 s. This marks the first instance that the 
evolution of the inner regions of an accretion disk has been observed in real-time. 

Subject headings: accretion, accretion disks — line: formation — stars: individ- 
ual (4U 1820-30) — stars: neutron — X-rays: binaries — X-rays: bursts 


1. Introduction 

Superbursts (e.g. Cornelisse et al. 2000; Strohmayer & Brown 2002; Kuulkers et al. 
2002) are the more powerful and rarer relations of the Type I X-ray bursts, and are thought 
to be due to the nuclear burning of the ashes (principally carbon and rp-process nuclei) 
that remain after thermonuclear processing of the accreted light elements (Strohmayer k. 
Brown 2002; Cumming 2003). One of the defining characteristics of a superburst is its long 
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decay time of a number of hours, three orders of magnitude greater than a typical Type I 
burst (Kuulkers 2003). This allows high signal-to-noise spectra to be accumulated at regular 
intervals throughout the course of the burst. 

4U 1820-30 is a compact low-mass X-ray binary with a period of only 11.4 minutes 
(Stella, Priedhorsky k White 1987), and lies within the globular cluster NGC 6624 at a 
distance of ~ 7 kpc (Vacca, Lewin k van Paradijs 1986; Hesser & Shawl 1985; Rich, Minniti 
&; Liebert 1993). Type I X-ray bursts have been observed from this system for nearly 30 
years (Grindlay k Gursky 1976). The accretion rate inferred from the persistent X-ray 
luminosity is ~ 1.5 x 10 17 g s _1 (Cumming 2003). The first, and so far only, superburst 
seen from 4U 1820-30 was observed by the Rossi X-ray Timing Explorer ( RXTE ) on 1999 
September 9 (UT) during a scheduled monitoring observation. The initial analysis and 
interpretation of the data was presented by Strohmayer & Brown (2002). Time-resolved 
spectroscopy presented by these authors showed that an iron Ka line and absorption edge 
were present in the spectra, and that their properties changed over time. A likely origin for 
these features is reprocessing of the X-rays by the accretion disk (Day k Done 1991). This 
process, called X-ray reflection, is commonly used to explain similar features in the spectra 
of Seyfert galaxies (e.g. Pounds et al. 1990; George k Fabian 1991; Ballantyne, Vaughan k 
Fabian 2003). 

In this paper we present the results of fitting the time-resolved spectra of the superburst 
with detailed models of X-ray reflection. The next section describes the model calculations 
and details of the spectral fitting. The results are presented in Section 3 and discussed in 
Section 4. 


2. Model Calculations 

The X-ray reflection models were calculated with the code of Ross k Fabian (1993) 
(see also Ross, Fabian k Young 1999; Ballantyne, Iwasawa k Fabian 2001; Ballantyne, 
Ross k Fabian 2001). A one-dimensional slab of uniform density gas is illuminated by 
a blackbody continuum (defined from 1 eV to 100 keV) of temperature kT, and allowed 
to come into thermal and ionization balance. The emergent spectrum is calculated with 
a Fokker-Planck/diffusion method that naturally accounts for the energy redistribution of 
photons by Compton scattering (Ross, Weaver k McCray 1978). The following ions are 
included in the calculations: C V - VII, N VI - VIII, O V - IX, Mg IX - XIII, S XI - XV 
and Fe XVI - XXVII. Fluorescence and recombination lines from these ions are computed 
and transferred simultaneously with the continuum, resulting in a self-consistent prediction 
of the emitted spectrum. 
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The compact orbit of this binary suggests that the donor is an evolved low-mass helium 
star (Rappaport et al. 1987), and thus the accretion disk material will be He-rich and H-poor 
(Cumming 2003). Therefore, we assumed the reflecting medium has the He, C, N, and O 
abundances of cool extreme helium stars from the tabulation of Pandey et al. (2001). The 
Mg, Si and Fe abundances were left at their solar values (Grevesse, Noels & Sauval 1996), 
but the effective metallicity is increased (by ~ 3) due to the negligible levels of hydrogen. 
The helium is assumed to be completely ionized everywhere in the reflecting region. 

The computed reflection spectra are parameterized by the ionization parameter £ = 
47rFx/nHe of the irradiated gas, where F x is the flux of the incident X-rays and nn e = 
10 15 cm -3 is the assumed number density of helium. Since the blackbody emission is also 
seen directly from the neutron star, the fitted spectrum is the sum of the blackbody and the 
reflected emission, weighted by a factor R, the reflection fraction. This parameter measures 
the relative importance of reflection in the spectrum, and can be related to the covering 
factor of the accretion disk. 

Strohmayer & Brown (2002) describes the RXTE data and the process of accumulating 
the ~ 80 64 s spectra during the burst. The preburst spectrum was used as the background, 
and was thus subtracted from the data during spectral fitting (Strohmayer & Brown 2002). 
This persistent emission (amounting to a luminosity of 1.8 x 10 37 erg s -1 ) is driven by the 
accretion onto the neutron star, so using this spectrum as a background implicitly assumes 
that the burst will not effect this underlying emission. As we see below, we cannot be sure 
that this does not happen; however, Strohmayer & Brown (2002) found little difference in 
their spectral fits between the preburst background and the detector background. Since the 
flux of the superburst varies from 20 x to 3x greater than the preburst emission, then the 
shape of the background will not greatly effect the fit results. 

The spectral fitting was performed with XSPECv.ll.2.0bp (Arnaud 1996). The model 
consisted of a grid of reflection spectra (read in using the atable command) and attenuation 
from neutral gas along the line of sight (using the wabs model), parameterized by the column 
density N^. The Galactic absorption toward 4U 1820-30 is 1.5 x 10 21 cm -2 (Dickey &; 
Lockman 1990). The reflection spectra were blurred by the effects of relativistic motion 
within the Schwarzschild metric (Fabian et al. 1989). The inner radius of emission (r- m ) was 
allowed to vary, but the emissivity was fixed as r -3 , where r is distance along the accretion 
disk, and the outer radius was frozen at 200 GM/c 2 , where M is the mass of the neutron 
star. The inclination angle of the disk to the line of sight was fixed at 30°. The fitted 
energy range was 3-40 keV, but was changed to 3-15 keV 7000 s into the burst because of 
background dominating at higher energies. The errorbars on the best fit parameters were 
computed using the 2 a uncertainties for one parameter of interest. 
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3. Spectral Fitting Results 

The results of fitting the reflection models to the superburst data are shown in Figure 1. 
We find that an ionized reflector provides a good description to the spectra for most of the 
burst. During the first 1000 s the spectra are well fit with a large £, indicating significant 
ionization of the accretion disk. Over the next ~ 500 s £ decreases by more than an order of 
magnitude. For times between ~ 2500 and 3000 s the spectra cannot be fit by the reflection 
model, except for the last three points, where only very loose constraints could be placed on 
£. Strohmayer & Brown (2002) showed the spectrum was hardening very rapidly during that 
period, which indicates that some rapid evolution of the system was ongoing, and is not well 
modeled by our equilibrium calculations. However, the recombination timescale for H-like 
Fe is trecomb ~ l/n e O!Fe(T). The slowest timescale occurs at the highest temperature, which 
is 10 7 K at the surface of the reflector, so that t r ecomb ~ 10 -4 s, assuming the recombination 
coefficient of Verner & Ferland (1996). This is the same order of magnitude as the disk 
dynamical time at 10 GM/c ?. The recombination timescale of He II is ~ 0.1 s, much longer 
than the dynamical time, but faster than the 64 s spectra examined here. Therefore, the 
assumptions of ionization equilibrium or fully stripped helium cannot be the cause of the 
poor fits in this region of the lightcurve. A more fundamental change in the X-ray source or 
reflecting geometry may be required. Finally, after the last Earth occultation passage, the 
system was revealed to have settled into a final decaying state, with £ ~ 100 erg cm s -1 , 
although one other non-equilibrium event took place at ~ 6000 s. 

While the disk inclination angle was fixed at 30° for the spectral fits, there is a previous 
constraint of 35°-50° found by Anderson et al. (1997) from modulations in the ultraviolet 
flux (due to X-ray heating of the secondary). Relativistic disk features are also dependent 
on the accretion disk inclination angle (Fabian et al. 1989), so can provide a consistency 
check on the reflection hypothesis. In Figure 2 we plot the 68%, 90% and 99% confidence 
contours for r; n and the disk inclination angle 5500 s into the superburst. This segment was 
chosen as it has a low reduced % 2 and strong reflection features. While the disk fit seems to 
prefer a lower inclination angle, it is consistent with the results of Anderson et al. (1997) at 
the 2 a level. 

The observed evolution of £ from highly to weakly ionized is to be expected since 
the observed flux decays over the course of the burst (Strohmayer & Brown 2002). More 
fundamental insights from this analysis arise from the variations of Ah, R and r; n over the 
first 2000 s of the superburst. The changes of both r; n , which measures the inner radius of 
the disk (as judged from the width of the iron line and edge), and R, the reflection fraction, 
strongly indicates that the accretion disk is significantly affected by the superburst out to a 
radius of ~ 100 GM/c 2 . In addition, the line-of-sight absorbing column also went through 
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a similar variation, either as a result of the burst itself or the changes to the inner accretion 
disk. All three of these parameters follow a similar pattern where they increase (r, n and 
Ah) or decrease ( R ) before returning to close to their initial value. This behavior seems to 
be closely linked to the blackbody temperature kT, rather than the flux, which decreases 
monotonically during the burst. 


4. Discussion 

A likely scenario to explain the observed behavior is that the accretion flow in the inner 
disk is altered by the superburst so that it loses the ability to produce reflection features. 
This could be accomplished if the scale-height of the disk, H, increased due to heating of the 
gas by the radiation field. This would decrease the surface density of the disk, allowing it to 
be more easily and completely ionized by the burst radiation, removing any spectral features 
in the observable band, and effectively becoming a mirror for the blackbody spectrum. As 
there would be no evidence for reprocessing from this innermost region, the estimate of the 
inner radius would be from further out in the disk where iron can still produce an emission 
line. The additional blackbody emission reflected from the ionized mirror, would dilute the 
reflection features, resulting in the low value of R. Alternatively, the lower surface density 
could reduce the optical depth of the gas to the point where the inner disk becomes a very 
inefficient reflector, while only the outer disk would be optically thick enough for efficient 
reprocessing. 

The main evidence for this thermal effect is seen in the correlations between r- m , R and 
kT. The Compton temperature of the gas is kTc ~ kT, and reached ~ 3 keV, about twice the 
maximum disk temperature (Shakura & Sunyaev 1973). Since H oc Cgr 3 / 2 oc T 1//2 r 3 / 2 , where 
Cg is the local sound speed, then H will increase by ~ 40% from its non-illuminated value. 
Changes to the disk surface density can take place over a viscous time which is ~ 1000 s 
at 100 GM/(? if a ~ 0.2, where a is the Shakura-Sunyaev viscosity parameter (Shakura & 
Sunyaev 1973). This timescale compares favorably with the observed changes in R and r- m . 
An increase in H may also explain the evolution of Ah as the thicker disk may bring clumps 
of disk material into the line-of-sight. Data which extend below the 3 keV limit of RXTE 
would be needed to test if the absorption was caused by ionized or neutral material. Under 
this hot disk scenario, once the peak in kT is passed the disk material cooled to lower and 
lower temperatures, and thus contracted to its original size. Simultaneously, the amount of 
material available to reflect increased, so R grew larger, but f dropped due to a combination 
of the increasing surface density and the decaying burst flux. 

The connection between the ionization parameter and the disk surface density E can be 
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formulated as £ « Lxmn e H/r 2 T,, where Lx is the luminosity of the burst, and nine is the 
mass of a helium nucleus. Using the equations for a gas pressure dominated Shakura-Sunyaev 
disk with electron scattering opacity (Shakura & Sunyaev 1973), a peak burst luminosity of 
3.4 x 10 38 erg s -1 (Strohmayer &; Brown 2002), and a = 0.2, we find £ « 100 erg cm s -1 for r = 
15 GM/c 2 , which is consistent with the values of £ found near the end of the burst (Fig. 1). 
From the observed decay in flux (Strohmayer & Brown 2002), the surface density of the disk 
must be ~ 100 times less than the Shakura-Sunyaev value at the start of the superburst, but 
only ~ 6 times less near the end in order to obtain the necessary level of ionization. This 
indicates the surface density of the inner accretion disk was significantly lower at the peak 
of the burst but increased as the burst decayed, supporting the interpretation that the disk 
material was severely affected by the explosion. 

Another process that could lower the surface density at the inner edge of an accretion 
disk is the interaction between the radiation field and material flowing in the disk. Radiation 
can remove angular momentum from the disk material, causing the disk to empty inward, 
thus lowering the surface density. There exists only a few calculations of the interaction 
between an accretion disk and the radiation field from a luminous central source (Walker 
& Meszaros 1989; Walker 1992; Miller & Lamb 1993, 1996). These models predict that the 
radiation will exert a significant torque on the accretion flow, which increases the accretion 
rate close to the central object. This effect may be enhanced in the case of 4U 1820-30 
due to the power of the superburst (as compared to the typical Type I burst considered in 
these models), as well as an interaction between the disk and the expanding shell that was 
ejected from the neutron star just prior to the superburst (Strohmayer & Brown 2002) - the 
radius-expansion phase (Lewin, Vacca &; Basinska 1984). It is difficult to rigorously compare 
this idea to the observations because detailed modeling of the disk-radiation interaction for 
a superburst or a radius-expansion burst has not been performed. 

By the end of the observation, the accretion disk has recovered from its disruption by 
the burst, and is reflective down to a radius of 10-20 GM/c 2 , which is curiously larger 
than the innermost stable circular orbit of 6 GM/c 2 . Since this is not an X-ray pulsar, the 
disk would not be truncated by the star’s magnetosphere. Rather, this may be evidence 
for the disk-radiation interaction described above, and the material within 10-20 GM/c 2 is 
no longer dense enough to provide a reflection signature. At these late times, the fits also 
require strong reflection, with R > 1, indicating that the disk covers a large solid angle 
from the X-ray source, and thus is thin relative to the radius of the neutron star or has 
a non-negligible flare at large radii. The latter possibility is supported by the reflection 
from r > 100 GM/c 2 observed early in the burst. If the inner disk expanded from the X-ray 
heating, then a significant disk flare might be necessary in order to allow the distant material 
to be irradiated by the blackbody emission. 
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Dipping behavior is seen in the X-ray lightcurve of 4U 1820-30 between 6000-7500 s into 
the burst (see Fig. 7 in Strohmayer k Brown 2002). This kind of variability in low-mass X- 
ray binaries is often explained as variable absorption from a two-component ionized medium 
close to the neutron star (Frank, King & Lasota 1987). The dipping observed from 4U 1820- 
30 could then be due to the reformation of the accretion disk, however, the timescale is 
incorrect, with the dipping seen ~ 2000 s after the probable recollapse of the disk. However, 
there is an interesting correspondence between the times when dipping is observed and an 
increase in iV H derived from spectral fitting (Fig. 1). This may indicate that the dipping is 
caused by partially ionized material that is slowly recombining, but farther away from the 
accretion disk. An obvious source for such material is the gas that was cast off from the 
neutron star during the radius expansion phase. This material could be at a large distance 
from the system (10 14 cm if ejected at the escape velocity). 


We thank E. Agol, L. Bildsten, E. Quataert, N. Murray k C. Matzner for helpful 
discussions, and A. Fabian for providing the relativistic blurring functions. DRB is supported 
by the Natural Sciences and Engineering Research Council of Canada. 
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Fig. 1. — Results of fitting the time-resolved superburst spectra of 4U 1820-30 with ionized reflection 
models. The top panel shows the evolution of the reflection parameters: the ionization parameter £, the 
reflection fraction, and inner radius of the disk, ri n . The lower panel describes the evolution of the absorption 
column, kT of the incident blackbody, and the reduced x 2 of the spectral fits. The reflection fraction, r in , 
and iVn seem to be closely related to the evolution of kT. Data were fit from 3-40 keV for the points to left 
of the dashed line, and from 3-15 keV for the points to the right. Time equals zero corresponds to about 
100 s after the start of the burst. The hatched region denotes the interval where the spectra were hardening 
rapidly, and were not well described by the equilibrium reflection model. 
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Fig. 2. — Joint 68%, 90% and 99% confidence contours on r in and the disk inclination 
angle for a spectral fit 5500 s into the superburst. Modulations in the ultraviolet flux from 
4U 1820-30 previously constrained the inclination angle to between 35° and 50° (Anderson 
et al. 1997), which is consistent with the results from disk reflection at the 2 o level. 




